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Continuous Transition between Antiferromagnetic Insulator and Paramagnetic Metal 

in the Pyrochlore Iridate Eu2lr207 
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Our single crystal study of the magneto-thermal and transport properties of the pyrochlore iridate 
Eu2lr207 reveals a continuous phase transition from a paramagnetic metal to an antiferromagnetic 
insulator for a sample with stoichiometry within ~1 % resolution. The insulating phase has strong 
proximity to an antiferromagnetic semimetal, which is stabilized by several % level of the off- 
stoichiometry. Our observations suggest that in addition to electronic correlation and spin-orbit 
coupling the magnetic order is essential for opening the charge gap. 



The effects of spin-orbit (SO) coupling have been one of 
the recent central subjects in condensed matter physics. 
Along this stream, extensive theoretical and experimen- 
tal work has been performed on topological insulators 
and elucidated the gaplcss helical surface state in weakly 
interacting semiconductors^'^. Significant roles of the 
SO coupling have been also found in various electron 
correlation effects in the 5d transition metal based ox- 
ides. In particular, the study on the tetragonal iridium 
oxide Sr2lr04'^'^ has indicated that the strong SO cou- 
pling splits the otherwise wide 5d t2g band into two bands 
based on an effective pscudospin joff = 1/2 doublet and 
a Jeff = 3/2 quadruplet. It has been proposed that the 
half-filled relatively narrow joff = 1/2 band at the Fermi 
level (^^f) opens the Mott gap due to electron correlation. 

For the study of the interplay between electron correla- 
tion and SO coupling in 5d transition metal compounds, 
pyrochlore oxides are suitable because the corner shar- 
ing tetrahedra network tends to form a narrow flat band 
at Ep that enhances the effects of electron correlation 
and/or SO coupling, both of which have comparable size 
with the bandwidth. In particular, rare-earth pyrochlore 
iridates i?2lr207 {R ~ rare-earth elements) provide an 
ideal class of materials because they have been reported 
to show a variety of transport and magnetic behaviors^"^, 
varying from a correlated metal^"^^ to an insulator^^'^^ 
depending on the ionic radius of rare-earth elements. 
They are also attractive because the one hole nature of 
Ir^+ 5d ilg state simplifies the electronic structure and 
allows further in depth comparison between experiment 
and theory. Indeed, a number of interesting theoretical 
proposals have been made as possible ground states of 
the iridium based pyrochlore oxides, such as a topolog- 
ical insulator, a Weyl semimetal, and an antiferromag- 
netic insulator^''"^®. 

To clarify the basic magnetic and transport proper- 
ties due to the iridium pyrochlore network, the simplest 
to study is the systems with rare-earth elements with 
nonmagnetic properties, namely, Eu^+ or Y'^^. While 
Y2lr207 is known to be insulating at all the temperature 
measured, the experimental situation for Eu analog is 
controversial. Pioneer work by Yanagishima and Maeno 
reports that it is a nonmagnetic metal^, while recent 
work by Matsuhira et al. found a metal-insulator transi- 



tion with a magnetic ordering in the insulating phase^'®. 
Recent Raman study found a subtle structural symme- 
try change^^ while the single crystal X-ray indicates no 
sign of the structural transition^". This suggests that 
Eu2lr207 has proximity to a metal-insulator transition 
and therefore best suited for the study of the interplay be- 
tween the electronic correlation and spin-orbit coupling. 

In order to reveal the intrinsic behavior of this interest- 
ing pyrochlore iridate, we have performed detailed single 
crystal study of transport and magnetic properties and 
their systematic dependence on sample quality. It reveals 
that the ground state of the stoichiometric Eu2lr207 is 
an antiferromagnetic (AF) insulator with a small charge 
gap (~ 10 meV) and with magnetic isotropy. The re- 
sistivity is exponentially divergent on cooling, exclud- 
ing the possibility of the topological semi-metallic state 
predicted by theory^^'^^ . Instead, this insulator has 
strong proximity to an AF semimetal or metal phase 
which may be stabilized by carrier doping induced by 
a few % off-stoichiometry. The fact that the onset of the 
AF ordering is independent of the magnitude of the off- 
stoichiometry indicates that the AF order is essential for 
opening the charge gap. These results are fully consistent 
with the recent theoretical proposal of the antiferromag- 
netic insulating ground state with "all-in all-out" spin 
structurc^^'^®. 

Single crystals of Eu2lr207 were grown by a KF flux 
method from polycrystallinc Eu2lr207, prepared by solid 
state reaction of the appropriate mixture of EU2O3 (4N) 
and Ir02 (4N) powders^". Powder and single crystal X- 
ray diffraction analyses confirm single phase and the py- 
rochlore structure with the lattice constant a = 10.274(3) 
A. Electron probe microanalysis (EPMA) found a slight 
deviation up to several % from the stoichiometry in the 
ratio between Eu and Ir contents. Our detailed study 
using synchrotron X-ray measurements of the related 
rare-earth iridate Pr2lr207 found the excess rare-earth 
/ iridium occupies the site of the counterpart ion with- 
out leaving vacancy^^. Therefore, throughout this paper 
we adopt the chemical formula Eu2(i_a;)Ir2(i-|_:i-)07+d-. As 
we will discuss below, our measurements indicate that 5 
should be zero, when x = Q. The change in the lat- 
tice parameter due to the off-stoichiometry was not de- 
tected within experimental resolution using our labora- 



tory based X-ray systems. This suggests that the major 
effect would be the carrier doping due to the valence im- 
balance between Eu'^+ and &■*+ ions. 

The longitudinal and Hall resistivity was measured by 
a standard four-probe method from 2 K to 300 K us- 
ing a variable temperature insert system under a field 
along [111] up to 9 T. DC magnetization under a field 
along [100], [110], and [111] axes was measured using 
a commercial SQUID magnetometer (MPMS, Quantum 
Design) from 2 K to 350 K under a magnetic field of 0.1 
T in both field cooled (FC) and zero field cooled (ZFC) 
conditions. Specific heat was measured by thermal relax- 
ation method using a commercial system (PPMS, Quan- 
tum Design). 

First we present the temperature dependence of the 
electrical resistivity p{T) . To discuss the strong sensitiv- 
ity to the off-stoichiometry. Figure 1 shows p{T) normal- 
ized by its room temperature value, p(T)/p(300 K), for 
single crystals with various off-stoichiometry x. In the 
high temperature region of T > 120 K, all the data of 
p(T)/p(300 K) collapse on top of each other and shows 
metallic behavior with positive slope with a typical value 
of 10 milcm at room temperature. The loffe-Regel limit 
was estimated to be ~ 1.4 mJ7cm, indicating that trans- 
port is strongly incoherent in this high T regime. 

On cooling in the range of 120 K > T > 100 K, all the 
samples commonly show a broad upturn in p(T) with an 
anomaly in dp/ AT at T = 120 K (inset of Fig. 1). Here 
at T = 120 K we locate the onset of the low temperature 




80 100 120 140 160 
r(K) 



100 150 200 

r(K) 



300 



FIG. 1. (Color online) Temperature dependence of the electri- 
cal resistivity of single crystalline Eu2lr207 for three samples 
with different off-stoichiometry x, normalized by the value at 
300 K. Inset: Temperature dependence of the temperature 
derivative of the resistivity dp/ AT normalized by the resistiv- 
ity at 300 K. 
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FIG. 2. (Color online) Inverse resistivity ratio 1/RRR 

= p(2 K)/p(300 K) vs. the off-stoichiometry |a::| for various 
single crystals of Eu2(i_i)Ir2(i+i)07+i. Blue and red symbols 
correspond to Ir rich (x > 0) and Eu rich (x < 0) samples, 
respectively. The blue shaded belt is a guide to the eye. 



nonnietallic (negative Ap/AT) behavior and refer to this 
as the metal- nonmetal transition temperature Tmn. No 
hysteresis in p(T) was observed across TmN; indicating 
the second-order nature of the transition. 

Below Tmn, piT) is found substantially different with 
various residual resistivity {p(T = 2 K)) varying by four 
orders of magnitude between various samples. To clarify 
the origin of the sample variation in p{T), we plot in Fig. 
2 the inverse resistivity ratio 1/RRR = p{2K)/p{300K) 
vs. the off-stoichiometry x for various single crystals. 
The deviation from stoichiometry is found at most 5 
% and the majority of samples are Ir rich. There is a 
correlation between \x\ and 1/RRR, which suggests that 
the most stoichiometric samples have the largest 1/RRR 
i.e. being the most insulating. This also indicates when 
\x\ — !• 0, S should also vanish, so that no doped carri- 
ers remain. Therefore, we conclude that the intrinsic low 
temperature state of Eu2lr207 is insulating and 1/RRR 
can be used as figure of merit for the sample quality. 

The insulating behavior can be also confirmed in the 
temperature dependence of the resistivity p(T). Figure 1 
shows that the most stoichiometric crystal =ffl with the 
smallest value of \x\ = 0.015(8) exhibits the most resis- 
tive behavior and a clear shoulder at ~ 50 K. As for the 
crystal #2 with larger x = 0.026(9), p{T) saturates at 
low temperatures, resulting in a small value of 1/RRR 
= p(2 K)/p(300 K). The sample with a smaller 1/RRR 
has a weaker anomaly in Ap/AT at Tmn, thus indicat- 
ing that off-stoichiometry broadens the transition. For 
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FIG. 3. (Color online) \np(T)/p(300 K) vs. l/T^^^ of sin- 
gle crystalline Eu2lr207 for three samples with different off- 
stoichiometry x. The solid lines are guides to the eye. Inset: 
Temperature dependence of the activation energy E^. For 
definition, see text. 



the sample with the largest x ~ 0.037(8) among three 
crystals used in Fig. 1, p{T) shows only a weak kink at 
Tmn and display a semi-metallic behavior with a gradual 
increase of the resistivity on cooling. To further char- 
acterize the insulating behavior, we plotted hip(T) vs. 
1/T" with various a. As shown in Fig. 3, the vari- 
able range hopping (VRH) p{T) = po exp(To/T)i/4 with 
a — 1/4: for three dimensions is found to best describe 
the insulating behavior^^"^^ , in particular, for the crystal 
#1 with a;=0. 015(8) over a decade of T below the shoul- 
der temperature of 50 K and down to the lowest T of the 
measurement, 2 K. With increasing the off-stoichiometry, 
p{T) deviates from the VRH fit (solid line in Fig. 3) and 
levels off at low temperatures. This strongly suggests 
that the insulating state of the near stoichiometric sam- 
ple with x=0. 015(8) has strongly localized states at the 
Fermi level Ep induced by the slight off-stoichiometry. A 
rough estimate of the charge gap can be made using the 
Arrhenius's law p{T) = po exp{Ea/kBT) with po being 
the minimum of the resistivity. The resultant T depen- 
dent Ea{T) for the crystal #1, shown in the inset of Fig. 
3, reaches its maximum of 16 meV at 40 K, indicating 
that Eu2lr207 is a narrow gap semiconductor. 

Comparing our single crystal results to recent results 
for polycrystalline samples* finds good consistency in the 
value of 1/RRR ~ 10^, providing another evidence for 
the high quality of our stoichiometric single crystal. More 
than ten times larger absolute value of the resistivity and 
the nonmetallic behavior with negative slope found even 



at T > Tmn in the polycrystalline samples indicate ex- 
trinsic scattering due to the polycrystalline grain bound- 
aries, which is absent in our single crystals. 

To further reveal the nature of the transition at Tmn, 
we measured the magnetoresistance and Hall resistivity 
as a function of field and temperature. Except a slight 
broadening of the kink at Tmn, no major change was ob- 
served in p{T) under field up to 9 T applied along [111], 
revealing the robust feature of the transition against the 
application of the field. The field dependence of the Hall 
resistivity p^y plotted in Fig. 4 shows a linear depen- 
dence on magnetic field poH up to at least 9 T in the 
temperature region between 90 and 140 K across Tmn- 
This linear response allows us to define the Hall coeffi- 
cient defined by Rn = pxy/noH, whose temperature de- 
pendence is plotted in the inset of Fig. 4. Interestingly, 
no anomaly was found at Tmn in sharp contrast with the 
longitudinal resistivity results. This behavior is similar 
to the one observed for Cd2 0s207^^. As we found above, 
the charge gap scale is only of the order of 10 meV com- 
parable with the scale of Tmn- Thus, thermal excitations 
across the gap should be significant at ~ Tmn, and may 
effectively wipe out the charge gap, smoothening the T 
dependence of Ra . 

Note that pxy at 90 K shows small hysteresis suggestive 
of time reversal symmetry breaking at T. Actually, hys- 
teresis in M-H curve (not shown) also appears just below 
Tmn with evolution of a small spontaneous magnetiza- 
tion. Hysteresis in pxy appears to come from the anoma- 
lous component proportional to the magnetization due 
to spin-orbit coupling and/or impurity scattering^^^^^. 
Considering the hysteresis in the MH curve together 
with the one observed for the susceptibility as a function 
of temperature described below, weak spontaneous mag- 
netization may originate from canting of non-collinear Ir 
moments. 

Figure 5 shows the temperature dependence of the 
magnetic susceptibility x{T) — M{T)/ H for a batch of 
5 single crystals with p(4.2 K)/p(300 K) > 500. The 
measurement was made using zero field cool (ZFC) and 
field cool (FC) sequences under a field of 0.1 T aligned 
along various high symmetry crystal axes. A sharp tran- 
sition was observed as a kink in x(T) at T = 120 K, the 
same temperature as Tmn- In addition, no hysteresis was 
observed across Tmn, indicating the transition is of the 
second-order type. The measurements for other batches 
of crystals with respective /o(4.2 K)/p(300 K) > 100 and 
with /9(4.2 K)/p(300 K) < 33 found the same onset tem- 
perature of 120 K, and in particular for the FC sequence, 
the nearly identical temperature dependence of xC^)- 
This is striking given the strong sensitivity of p{T) to 
the off-stoichiometry. The magnetic transition at Tmn is 
also in good agreement with the observation of a mag- 
netic order observed by muon spin rotation (pSR)^^. 

Above Tmn, xC^) shows a weak temperature depen- 
dence. Generally for Eu compounds, it is known that the 
Van-Vleck paramagnetism of Eu^+ ions has to be consid- 
ered because there are relatively low lying J — \ {^ 480 
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FIG. 4. (Color online) Field dependence of the Hall resis- 
tivity pxy for single crystalline Eu2lr2 07 with 1/RRR = 200 
obtained under external field along [111] at various tempera- 
tures between 90 K and 140 K across the transition tempera- 
ture of Tmn ~ 120 K. Inset: Temperature dependence of the 
Hall coefficient Rh = pxy/B obtained under a field of 9 T 
applied along [111]. 



K) and J = 2 (-- 1440 K) excited multiplets^^ . Figure 
5 shows the Van-Vleek eontribution xvv estimated us- 
ing the above exeited multiplet levels for a J = ground 
state. Note that the contribution obtained after subtract- 
ing the Van-Vleck contribution from the observed x{T) 
is nearly temperature independent and does not follow 
the Curie- Weiss form. This strongly suggests that this 
component comes from the Pauli paramagnetism and in- 
dicates the absence of local moment in the metallic phase. 

At Tmn = 120 K the ZFC data exhibit a sharp 
peak followed by a drop in x(^) on cooling. This 
drop exists for crystals with p(4.2 K)/p(300 K) > 500 
and /9(4.2 K)//9(300 K) > 100, but not for those with 
p(4.2 K)/p(300 K) < 33 or the results reported for poly- 
crystalline samples^, wherein x{T) returns to the value 
for the metallic phase. Because the size of xC^) after 
the drop for the high 1/RRR samples is close to the one 
for Xvv, a possible explanation for this drop is the dis- 
appearance of the Pauli paramagnetic component as the 
system becomes insulating by opening a gap at E-p- We 
note the temperature range in which this drop occurs is 
concurrent with the initial increase in p at T < Tmn- 
The absence of the drop suggests no well defined charge 
gap, as is consistent with the low value of 1/RRR. 

For both FC and ZFC results, the magnetization is al- 
most independent of the magnetic field alignment, indi- 
cating that unit cell of the spin structure should be mag- 
netically isotropic. This suggests a non-coUincar arrange- 



FIG. 5. (Color online) Temperature dependence of the mag- 
netic susceptibility x ~ M/H of single crystalline Eu2lr2 07 
with 1/RRR > 100, measured using zero field cool (ZFC) and 
field cool (FC) sequences under a field of 0.1 T aligned along 
[100], [111], and [110]. The broken lines indicate the Van- 
Vleck contribution xw estimated using the excited multiplet 
levels for a J = ground state of Eu^"*" ion and the compo- 
nent obtained after subtracting the Van-Vleck contribution 
from the observed x{T), labeled as Xobs — Xvv. Inset: x{T) 
for three difi'erent qualities of crystals under a field of 0.1 T 
along [110]. Both FC and ZFC results are shown. 



ment of Ir moments, and in particular, together with the 
commensurate nature found by /^SR measurements'^^, the 
"all-in all-out" antiferromagnetic ground state predicted 
theoretically for the pyrochlore iridates^^'^^. 

At low temperatures below 5 K, the upturn is observed 
in x{T) for the samples with p(4.2 K)//9(300 K) > 500 
and p(4.2 K)/p(300 K) > 100, while a drop between 
3 and 5 K for the one with p(4.2 K)/p(300 K) < 33. 
Nearly the same magnitude of the upturn was also ob- 
served in the results obtained for our polycrystalline sam- 
ple and suggests that this comes, not from the surface, 
but from the magnetic impurity in the bulk induced by 
ofF-stoichiometry. 

Figure 6 shows the temperature dependence of the 
specific heat divided by temperature Cp/T for samples 
with p(4.2 K)/p(300 K) > 500. Cp/T exhibits a clear 
anomaly at Tmn = 120 K providing evidence of the 
bulk transition. The measurements for samples with 
/9(4.2 K)/p(300 K) < 33 did not detect any anomaly in 
Cp and suggests that the off-stoichiometry may broaden 
the transition, as seen in the T dependence in dp/dT 
(inset of Fig. 1). 

At low temperatures we estimate the electronic specific 
heat coefficient, 7, by plotting Cp/T vs. T^. A linear 
fit is applied between 3 and 6.3 K yielding 7 — 12.9(1) 
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FIG. 6. (Color online) Temperature dependence of the spe- 
cific heat divided by temperature, Cp/T, of the single crys- 
talline Eu2lr207 with 1/RRR > 500. Inset: Cp/T vs. T^ for 
crystals with 1/RRR > 500 and with 1/RRR < 33. The solid 
lines are guides to the eye. 



mJ/niol-K^ and 7 = 24.9(2) mJ/mol-K^ for samples with 
/9(4.2 K)/p(300 K) > 500 and with p(4.2 K)/p(300 K) < 
33, respectively. The slope of the fits are the same and 
gives the Debye T to be 472(2) K. As expected, more 
conductive samples have larger 7 due to the presence 
of additional carriers. The finite 7 for insulating sam- 
ples could be associated to the strongly localized states 
at E-p induced by residual disorder. By further assum- 
ing a spherical Fermi surface for the localized states to 
estimate the density of states N{Ep) at Ep, the localiza- 
tion length ^ may be calculated^'^"^'^ , using the equation 
£_ = [21/(fcBroiV(-BF))]^/^- This yields the estimate of 
^ ~ 5 A, which is less than twice of the Ir-Ir ionic dis- 
tance, consistent with the strong localization. 

Below T = 5 K we find a striking upturn in Cp/T 
for both insulating and metallic samples. In the same 
temperature region, the upturn in x(T) was observed and 
therefore, this specific heat anomaly may be attributed 
to the impurity moments. In particular, for the sample 
with /3(4.2 K)/p(300 K) < 33, a kink is found at around 
3 K in Cp{T), similarly to the one found in x(T). This 
suggests a magnetic phase transition induced by a high 
concentration of the impurity spins. 

Summarizing all the results based on our single crystal 
study, we find the continuous, most likely second-order, 
phase transition at Tmn = 120 K from the paramag- 
netic metal to the antiferromagnetic insulator, which is 
intrinsic to the stoichiometric Eu2lr2 07. The gap size es- 
timated from the activation law is of the order of 10 mcV, 
which is the same size of the magnetic ordering tempera- 
ture, Tmn- This low temperature transport changes dra- 



matically from the insulating behavior to semi-metallic 
one, strongly depending on the sample stoichiometry. In 
sharp contrast, the magnetic transition temperature and 
the temperature dependence of the susceptibility is not 
affected by the off-stoichiometry. This clearly indicates 
the strong proximity of the ground state to the metal- 
insulator transition, and that the magnetic transition is 
essential for the system to fully open the charge gap when 
it is stoichiometric. 

According to the recent theoretical calculations for the 
pyrochlorc iridates, various types of interesting ground 
states have been predicted depending on the strength of 
the correlation and transfer integrals between the neigh- 
boring Ir 5d orbitals^^"^*. Among them, the phases 
with realistic values of these electronic parameters con- 
tain a metal, a topological semimetal, and an antiferro- 
magnetic insulator. Furthermore, both the calculations 
based on the local spin density approximation (LSDA) 
as well as the tight binding approximation predicts the 
magnetic structure to be the "all-in all-out" type struc- 
ture reflecting the strong spin-orbit coupling due to Ir 5d 
electrons^^'^*. Our observation of magnetic isotropy be- 
low Tmn as well as a commensurate order found by /xSR 
measurements'^^ are consistent with the "all-in all-out" 
magnetic structure. 

Strong excitement made by the theoretical predic- 
tion is actually for the possible existence of the topo- 
logical semimetal phase with Weyl or two-component 
fermions^^'^^. It has been proposed that the topologi- 
cally protected Dirac points are stabilized, which can be 
viewed as a three dimensional analog of the electronic 
structure of graphenc, and have an interesting conse- 
quence of a metallic surface state. ^^'^^ For this partic- 
ular semi-metallic state, the recent detailed calculation 
indicates that the resistivity should vary as a linear func- 
tion of 1/T and thus Tp{T) should be constant, reflecting 
the linearly vanishing density of states nearby the Fermi 
level. '^^ However, our results of p(T) is exponentially di- 
vergent on cooling, and excludes the possibility of such 
an exotic semi-metallic ground state. Instead, all our ob- 
servations are fully consistent with the prediction for the 
antiferromagnetic insulator with the "all-in all-out" mag- 
netic structure. Furthermore, the theory also suggests 
its proximity to a magnetic metallic phase with the same 
"all-in all-out" magnetic structure, which is also consis- 
tent with our observations. In addition, the small size of 
the gap is also the origin why the system is so sensitive 
to the off-stoichiometry and becomes (semi)metallic by 
the associated carrier doping. 

Finally, our recent study of the resistivity under pres- 
sure for Eu2lr207 found a quantum phase transition to a 
diffusive metallic phase from the AF insulator at the am- 
bient pressure. ^■^ According to the theoretical prediction, 
it is likely that the diffusive metallic phase under pres- 
sure may well correspond to the topological semi-metallic 
phase. Further study on the properties of this state under 
pressure is highly desired to elucidate the exotic nature 
of the ground state. 



We thank Yo Machida and Yasuo Ohta for the con- 
tribution in the early stage of this study. We also ac- 
knowledge L. Balents, M. Imada, S. Julian, Y. B. Kim, 
S. Onoda, and M. Takigawa for useful discussion. This 
work is partially supported by Grant-in-Aid for Scien- 
tific Research (No. 21684019) from JSPS, by Grant- 
in-Aid for Scientific Research on Priority Areas (No. 



19052003) from MEXT, Japan, by Global COE Programs 
"the Physical Sciences Frontier" MEXT, Japan, and by a 
Toray Science and Technology Grant. E. O.F. acknowl- 
edges the Japan Society for the Promotion of Science 
(JSPS) Fellowship for Foreign Researchers. The authors 
also thank the Materials Design and Characterization 
Laboratory, Institute for Solid State Physics, University 
of Tokyo. 



* jun@issp.u-tokyo.ac.jp 

t satoru@issp.u-tokyo.ac.jp 

^ M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 

(2010). 
^ X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 

(2011). 
3 B. J. Kim, H. Jin, S. J. Moon, J.-Y. Kim, B.-G. Park, 

C. S. Leem, J. Yu, T. W. Noh, C. Kim, S.-J. Oh, J.-H. 

Park, V. Durairaj, G. Cao, and E. Rotenberg, Phys. Rev. 

Lett. 101, 076402 (2008). 

* B. J. Kim, H. Ohsumi, T. Komesu, S. Sakai, T. Morita, 
H. Takagi, and T. Arima, Science 323, 1329 (2009). 

^ D. Yanagishima and Y. Maeno, J. Phys. Soc. Jpn. 70, 2880 
(2001). 

® N. Taira, M. Wakeshima, and Y. Hinatsu, J. Phys.: Con- 
dens. Matter 13, 5527 (2001). 

K. Matsuhira, M. Wakeshima, R. Nakanishi, T. Yamada, 
A. Nakamura, W. Kawano, S. Takagi, and Y. Hinatsu, J. 
Phys. Soc. Jpn. 76, 043706 (2007). 

* K. Matsuhira, M. Wakeshima, Y. Hinatsu, and S. Takagi, 
J. Phys. Soc. Jpn. 80, 094701 (2011). 

^ S. Nakatsuji, Y. Machida, Y. Maeno, T. Tayama, T. Sakak- 

ibara, J. vanDuijn, L. Balicas, J. N. Millican, R. T. 

Macaluso, and J. Y. Chan, Phys. Rev. Lett. 96, 087204 

(2006). 
^° Y. Machida, S. Nakatsuji, Y. Maeno, T. Tayama, T. Sakak- 

ibara, and S. Onoda, Phys. Rev. Lett. 98, 057203 (2007). 
^^ Y. Machida, S. Nakatsuji, S. Onoda, T. Tayama, and 

T. Sakakibara, Nature 463, 210 (2010). 
^^ M. Sakata, T. Kagayama, K. Shimizu, K. Matsuhira, 

S. Takagi, M. Wakeshima, and Y. Hinatsu, Phys. Rev. 

B 83, 041102 (2011). 
^■' F. F. Tafti, J. J. Ishikawa, A. McCollam, S. Nakatsuji, and 

S. R. Julian, arXiv:1107.2544vl [cond-mat.str-el]. 
" D. Pesin and L. Balents, Nat. Phys. 6, 376 (2010). 
^^ B.-J. Yang and Y. B. Kim, Phys. Rev. B 82, 085111 (2010). 
^^ X. Wan, A. M. Turner, A. Vishwanath, and S. Y. 

Savrasov, Phys. Rev. B 83, 205101 (2011). 
^"^ M. Kurita, Y. Yamaji, and M. Imada, J. Phys. Soc. Jpn. 



80, 044708 (2011). 
^* W. Witczak-Krempa and Y. B. Kim, Phys. Rev. B 85, 

045124 (2012). 
^' T. Hasegawa, N. Ogita, K. Matsuhira, S. Takagi, 

M. Wakeshima, Y. Hinatsu, and M. Udagawa, J. Phys.: 

Conf. Ser. 200, 012054 (2010). 
^° J. N. MiUican, R. T. Macaluso, S. Nakatsuji, Y. Machida, 

Y. Maeno, and J. Y. Chan, Mater. Res. Bull. 42, 928 

(2007). 
^ E. Nishibori, H. Sawa, Y. Ohta, and S. Nakatsuji, Preprint 

(2012). 
^^ S. Nakatsuji, V. Dobrosavljevic, D. Tanaskovic, M. Mi- 

nakata, H. Fukazawa, and Y. Maeno, Phys. Rev. Lett. 93, 

146401 (2004). 
■^■^ N. F. Mott and E. A. Davis, Electronic Processes in Non- 
Crystalline Materials (Clarendon Press, Oxfor, 1979). 
^* B. I. Shklovskii and A. L. Efros, Electronic Properties of 

Doped Semiconductors (Springer- Ver lag,, 1984). 
'^^ N. F. Mott, Conduction in Non-Crystalline Materials 

(Clarendon Press, Oxford, 1987). 
^^ D. Mandrus, J. R. Thompson, R. Gaal, L. Forro, J. C. 

Bryan, B. C. Chakoumakos, L. M. Woods, B. C. Sales, 

R. S. Fishman, and V. Keppens, Phys. Rev. B 63, 195104 

(2001). 
'^'^ N. Nagaosa, J. Sinova, S. Onoda, A. H. MacDonald, and 

N. P. Ong, Rev. Mod. Phys. 82, 1539 (2010). 
^* R. Karplus and J. M. Luttinger, Phys. Rev. 95, 1154 

(1954). 
^^ J. Smit, Physica 21, 877 (1955). 
^° J. Smit, Physica 24, 39 (1958). 
^1 L. Berger, Physica 30, 1141 (1964). 
^2 S. Zhao, J. M. Mackie, D. E. MacLaughlin, O. O. Bernal, 

J. J. Ishikawa, Y. Ohta, and S. Nakatsuji, Phys. Rev. B 

83, 180402 (2011). 
^^ Y. Takikawa, S. Ebisu, and S. Nagata, J. Phys. Chem. 

Solids 71, 1592 (2010). 
^^ P. Hosur, S. A. Parameswaran, and A. Vishwanath, Phys. 

Rev. Lett. 108, 046602 (2012). 



